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Hyper-fast interstellar travel... 


• Is there a way within the framework of physics such that one 
could cross any given cosmic distance in an arbitrarily short 
period of time, while never locally breaking the speed of light 
(11 th commandment)? 

SPACEWARPS 



Images courtesy NASA 



Inflation: Alcubierre Metric 1 



Warp Metric: 
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Alcubierre, M., "The warp drive: hyper-fast travel within general relativity," 
Class. Quant. Grav. 11, L73-L77 (1994). 
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Bubble Topology Optimization 

York Time magnitude decreases 





"bubble" thickness decreases 


Energy density magnitude decreases 



"bubble" thickness decreases 
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Surface plots of York Time & T°°, <v>=10c, 10 meter diameter volume, variable warp "bubble" thickness 

















































Bubble Oscillation Optimization 
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Exotic Mass Warp Requirements, 10m diameter, v_ QDarent =10c 
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Topology 


Approximations 
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Warp Field Interferometer 



• Warp Field Interferometer developed after 
putting metric into canonical form 1 : 

ds 2 = (y 2 f(r s ) 2 -1)< dt — > ~ dx 2 +dy 2 +dz 2 

. /W -! , 

• Generate microscopic warp bubble that 
perturbs optical index by 1 part in 
10,000,000 

• Induce relative phase shift between split 
beams that should be detectable. 



with device off with device on 
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White, H., "A Discussion on space-time metric engineering/' Gen. Rel. Grav. 35, 2025-2033 (2003). 
















Warp Bubble Detection Attempts Using 

Interferometry 


Goal is to use interferometer to detect and measure 
effect of warp bubble on optical path length through 
the measurement of associated interference fringe 
shifts. 

The EWL has attempted to mitigate the effect of 
vibrations & air currents 

— Using a vibration-isolated optical table 

— Using a vibration-isolated room 

— Using an optical hood 

— Using signal processing to increase signal to noise ratio 

— Collecting statistical data to increase signal to noise ratio 




Warp field 
Interferometer 





Interferometer and Test-article Setup 



2 Polarizers 

(jointly used for 

intensity 

control) 
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Primary/Auxiliary 



Isolator pressure 


Height control pressure 

















SINGLE PATCH 

INTERFEROMETER MODAL ANALYSIS, 
ANALYTIC SIGNAL VARIATION 
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Interferometer Modal Analysis: 

Synthetic Data _ 


1 


A set of 1000 synthetic interferogram fringe patterns were generated as JPEG image 
frames 

- Set contained repeating sequence of 10 frames with sinusoidal shifting in time with a 
shifting amplitude of 1/10 of a spatial radian*. 

- No noise added to the shifting (i.e. unwanted vibrations not included) 

- Gray scale intensity sampled at upper left pixel and a Discrete Fourier Transform (DFT) 
calculated for the sample train. 

- Very strong peak observed at 1/5 Nyquist frequency as expected, and of course, no noise 
spectrum. 

DFT Magnitude Spectrum 


Fringe pattern repeatedly 
shifted up and down along 
this axis 


\ 



Magnitude 


Strong peak observed at the 
excitation frequency (1/5 
Nyquist frequency) 


Excitation Frequency 


i 


Nyquist Frequency! 


2 fringe patterns out of set of 1000, showing maximum shift amplitude 
(In this case, the fringe pattern on the right is shifted about 1 radian up 
and to the left compared to the fringe pattern on the left.) 
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Modal Analysis: Analytic Signal Variation 

Single Patch 


2D DFT of the dashed region 



The result is an array of complex values, called the analytic signal. Now, set all DFT coefficients to 

Use the complex value at the center and calculate the "phase". zero except for this small 

Repeat for all of the captured frames and store the representative ^ - region around this peak 

phase values as a ID series. corresponding to the fringe 

lines, and inverse transform^ 
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Single Patch Modal Analysis, 

always off 



Magnitude Spectrum of Phase Series 
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Single Patch Modal Analysis, 

cycled on/off 



•Test article cycled on and off at .25 Hz, and associated 
fringe images frame sampled at 2.5 Hz for 3149 samples. 
•Any effect should show up at frequency 315. 

•Nothing around that frequency stands out above the noise 
for this run. 

•No evidence that the test article caused fringes to shift at 
this measurement precision. 

•Concern that single patch may not sample enough of 
image, so subsequent effort done to sample across image. 





Magnitude Spectrum of Phase Series 
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MULTIPLE SAMPLE POSITIONS 

INTERFEROMETER MODAL ANALYSIS, 
ANALYTIC SIGNAL VARIATION 


17 




Description of Work 


• Room and table were floating 

• Optical hood covering table 

• Test article cycled on every 4 seconds 

• Frames captured every 0.4 seconds 

• 3149 frames processed 

• Sampling occurred at 88 positions in image plane 

• At each pixel sampled (through 3149 frames), the 
sample train was Fourier transformed, and the 
magnitude spectrum reviewed (see next slides) 
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Magnitude 


Test Article Off 



Magnitude Spectrum of Phase Series 



Frequency 








Magnitude 


*c*o p«o**V 


Test Article Cycled every 4 
(2 sec @ 20KV, 2 sec off) 



Magnitude Spectrum of Phase Series 










Close up (off) 
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Close up (cycled) 
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Air Dielectric in Beam Path 


• Cap ring replaced with high-voltage device "Ozona" such that 
the charge collection sphere partially blocked the interferometer 
beam in the leg along the laser axis. 

• High voltage source was cycled on and off in 2 second intervals. 

• The screen capture was operating at a rate of 0.4 sec per frame. 

• The optics table was 
floating but not the 
room. 

• The optics table was 
not covered. 

• HV source power 
strip was on the 
micro flat table next 
to the interferometer 
table. 
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HV Power Supply Off 
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HV Power Supply Cycled @ 2 second 

Intervals 
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DAN NEHLICH, SOUTH DAKOTA STATE UNIVERSITY 

INTERFEROMETRY MEASUREMENTS 
OF A DC TOROIDAL CAPACITOR 
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SDSU Interferometer Setup 



Warp Field Interferometer 



DC Toroidal Capacitor Test Article 



Interference Pattern 
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f Image Subtraction Approach 

Positing that aberrations may include fringe pattern shifting in 
addition to or instead of induced curvature in the patterns, image 
subtraction was used to determine possible correlations to the 
voltage differential. 

- This method treats each image as an array and subtracts that array 
from the preceding, resulting in an image representative of the shift 
that occurred. 

- Average intensity of each subtracted image is calculated and statistical 
comparison is performed. 

- Statistical analysis of mean intensity values of each subtracted image, 
showed the observed shifts tended to increase during 
charge/discharge cycles, and remain constant in the control runs. 

The data presented in these scatter plots represents 2200 unique 
data points for charge/discharge cycles at 19 kV. 

Each point plotted represents approximately 20 images taken at the 
same point of the cycle, yielding an average intensity for each point 
in time along the graph. 



Control, Charge, Discharge Scatter Plots 



control 


charge 


discharge 


d(f) 

dt 
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Figure a: Image Subtraction! during Test Article Discharge Cycle from IS kl/to D V 
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Hough Transform Approach 


In order to establish a direction of fringe 
pattern shifting, edge detection 
algorithms were necessary to determine 
individual line locations. 


The Hough Transform, along with the 
associated houghpeaks and houghlines 
functions in Matlab, was utilized to this 
end. 


Hough Line Representation 


An algorithm was written to utilize the 
edge mapping function of the transform 
determine the distance of each line to 
the origin (normal to the fringe 
orientation) and then sort and weight 
the data to provide the relevant fringe 
locations for analysis. 


Hough Transform with Peaks 


Average Shifts of Interference Pattern 


Cycle 

Average Shift 
(pixels/frame) 

St. Dev. 

Average Shift -limited 
( pixel s/frame) 

St. Dev. 

Control 

0.050652543 

1.652173 

0.006819178 

0.319178 

Charge 

-0.012330713 

1.720073 

0.010439022 

0.341535 

Discharge 

0.005357319 

0.810996 

0.015817141 

0.188S45 














ENHANCED WARP FIELD PHYSICS EXPERIMENTS 

IMAGE AVERAGING, FABRY-PEROT 
INTERFEROMETER, & TIME OF FLIGHT 
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Fabry-Perot Warp Field Interferometer 


Front 

Mirror 


Imager 


DC or (preferably) AC 
Test Article 


Consists of two reflecting, 
highly parallel surfaces, called 
an Etalon 

The interference pattern is 
created within the Etalon 

Multiple reflections in the 
Etalon reinforce the areas 
where constructive and 
destructive interference 
occurs 

Allows for much higher- 
precision measurements of 
fringes (image averaging 
without software) 


Example: Michelson-Morley 
Interferometer image for 
Sodium source 


Example: Fabry-Perot 
Interferometer image for 
Sodium source (note doublet) 











SA-210 Scanning Fabry-Perot 
Interferometer 





Time of Flight Schematic 



200MHz test device 



Detectors Thorlabs High-Speed Avalanche Detector 



Model 

APD210 

Rise Time 

0.5 ns 

Supply Voltage 

+12 to +15 V 

Current Consumption 

200 mA 

Max. Incident Power 

10 mW 

Spectral Range 

400-1000 nm 

Frequency Range 

1-1600 mHz 

Maximum Gain 

2.5xl0 5 V/W 



HeNe 


EO Modulator 


Newport Broadband Amplitude Modulator 

Model 4102NF 

Type Broadband Amplitude Modulator 


Operating Frequency DC-200 MHz 


Wavelength Range 
Material 
Maximum Vk 


500-900 nm 
MgO:LiNbO s 
195 V <§> 633 nm 


Maximum Input Power 2 W/mm 2 @ 532 nm 


Aperture Diameter 

2 mm 

RF Bandwidth 

200 MHz 

RF Connector 

SMA 

Input Impedance 

10 pF 

Maximum RF Power 

10W 

Connector 

SMA 




Agilent Technologies Infiniium DS09254A 2.5 GHz Oscilloscope 

•2.5 GHz bandwidth across all 4 analog channels 
•20 GSa/s max. sample rate 

•Standard 20 Mpts memory per channel, upgradeable to 1 Gpts 
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Forward Plan 


DC Toroidal Capacitor Approach: 

- Work with larger sample sets to decrease effects of vibrational "noise" and 
develop test articles with longer regions of optical influence to increase the 
signal. 

- Utilize image averaging algorithms to increase fidelity of interference 
information. 

- Utilize the Fabry-Perot Interferometer to increase the sensitivity of the 
experimental apparatus to below 1/100 of a wavelength. 


Explore the d<p/dt dependency in future test devices 

- The idea of an optimized space warp needs negative vacuum energy, and large 
d<f)/dt - both of these conditions are present in the q-thruster technology also 
being explored in the lab. 

- Use the q-thruster physics models to guide design of RF frequency test devices 
to be evaluated in the warp field interferometer, the Fabry-Perot 
Interferometer, and the time of flight experiment. 




Casimir force 


Q-th rusters are a low-TRL form of electric propulsion that operates on the 
principle of poshing off of the quant oro vacu u m. 

A terrestrial analog to this Is to consldler how a submarine uses its propeIler to 
posh a colum n of water in one direction* while the sob recoils in the other to 
conserve momentum - the submarine does not carry a "tan k" of sea water to be 
used as propellant. 




Quantum 

Vacuum 

MHD 


In oor case* we ose the tools of Magnetohydrodynamics (MHD) to show how the 
thruster poshes off of the quantum vacuum which can be thought of as a sea of 


existence* and where fields are stronger* there are more virtual particles. 

The idea of pushing off the quantum vacuum has been in the technical literature 
for a few decades* but to date* the obstacle has been the magnitude of the 
predicted thrust which has been derived analytically to be very small* and 
therefore not likely to be useful for human spaceflight, 

©or recent theoretical model development and test data suggests that we can 
greatly increase the magnitude of the negative pressu re of the qua ntum vacuum 
and generate a specific force such that technology based on this approach can be 
competitive for In-space propulsion ('“O.lN/kW)* and possibly for terrestrial 
applications (~10N/kW). 

As an additional validation of the approadto * the theory a 11 o ws calculation of 
physics constants from first principles: Gravitational constant* Planck constant* 
Bohr radius* dark energy fraction* electron mass. 











Q-thruster Physics Data 


2004 Test Article 


£JJ * » ifr 


2005 Test Campaign 


2012 Test Article 




The test unit was run at 2.13 MHz, yielding an AC electric field 
of “20kV/m, and an AC magnetic field of “27 Gauss. 


• The observed thrust was +/- 0.89 mN 


The test unit was run at 3.8 MHz yielding an AC electric field 
of “20kV/m, and an AC magnetic field of “48 Gauss. 



• The observed thrust was +4.91 tc 
4900 mN (500gf) load cell 


-1.96 mN as measured vi 


As can be seen to the right, the thrust signal is very clear 
when the unit is excited. 


i uipws 


2012 test article tested in November 
2012 

2012 test article tested -> 98uN 
predicted, 2-3 uN detected 

• Scientifically very significant as vacuum 
fluctuation density had to be increased from 
-1x10-26 to > 1x10-14 

• As built quality factor much lower than 
desired, more engineering work necessary 

• Adjustment to power distribution network 
are in work to address power losses, 
increase thrust level 



SFE Test Article at JSC 


In 2013, Boeing/DARPA sent Eagleworks Lab an SFE 
test article for testing and evaluation 

Evaluation of the test article in and out of a Faraday 
Shield performed from Feb through June 2013. 

• There is a consistent transient thrust at device 
turn-on and turn-off that is consistent with Q- 
thruster physics 

• The magnitude of the thrust scaled 
approximately with the cube of the input 
voltage (2Q-110uN}. 

• The magnitude of the thrust is dependent on 
the AC content of the turn-on and turn-off pulse 

• Specific force of transient thrust was in the ~1- 
20 N/kW range. 

-20-110 uN Thrust Pulses 
Specific Force -1-20N/kW 



-3 mN Thrust 
Specific Force ~0.3N/kW 


Microwave Thruster Device 


SPR Ltd. Has produced several Microwave test articles. Claim is they produce 
thrust with just electric power input. . 

• Shawyer's theoretical model has been deemed 
non-viable by scientific community ( rightly so l 

Thruster assessed against Q-thruster models and 
w.” analysis suggests this may be a microwave 

version of a quantum vacuum plasma thruster. 

• Tapered shape creates virtual toroid of active volume Prototype 16mN @ 850W, 

that can realize net thrust in virtual plasma. 0.02N/kW 

• Microwave Q-thrusters would not be restricted to 
tapered construction. 

-4 • Thrust magnitude increased over multiple test 
devices from 16 to 170mN 
If Q-Thruster theory accounts for measured 
force, then microwave test articles may 
have ability to reach >10N/kW 
Chinese university claims to have duplicated 
EM Drive tests, but no way for U.S. to 
evaluate credibility (so we have Ignored itj 


[ 


-vacuum fluctuation density 
increased from 1x10 _26 to >1x10‘ 




* 




Cannae Test Article 



At DARPA's request, Eagleworks Lab began discussions with 
Cannae LLC on hosting a test article 
Cannae's "Q-Drive" design is inspired by RF resonant cavity 
design used in particle linear accelerators 

• Cannae theoretical model is likely non-physical. 

Initial assessment of conditions in thruster during operation 
indicate it may be a high frequency q-thruster. 


16-170 mN Thruster 
Specific Force 0.02-0.4N/kW 




High fidelity Test Article 170mN 
@ 450W, 0.4N/kW 


-7-10 mN Thruster 
Specific Force ~0.75N/kW 



Second Generation Test Article 
(under development) 


What we suspect: 

• A variety of industry experiments, for which theory is lacking, may be Q-thrusters 
including Boeing, Lockheed-Martin, EM Drive, Cannae, etc. 

• Low measured thrust but specific power ranges from 0.3 to 10+ N/KW 
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Q-thruster Roadmap 
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Q-THRUSTER + 2MW NUCLEAR POWER IS MISSION ENABLING 

POSSIBLE MISSIONS TO MARS, THE 
OUTER SOLAR SYSTEM, AND BEYOND 
WITH Q-THRUSTERS 
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Mars 




90 day stay 



• 140 day mission with 90 day 

stay at Mars 
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Jupiter 



• 2MW power 
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ARTH 


0.4 N/kW 
263 days 


SATURN 


SATURN 


4 N/kW 
86 days 


• 

90t spacecraft 

J^^ARTH 

• 

- 50t cargo, 20t power, 
20t propulsion 

2MW power 
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0.4 N/kW 
399 days 


JJ RAMUS 



4 N/kW 
129 days 


90t spacecraft 

- 50t cargo, 20t power, 
20t propulsion 


• 2MW power 
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Neptune 




APJH • 90t spacecraft 

- 50t cargo, 20t power, 
20t propulsion 

• 2MW power 


NEPTUNE 


NEPTUNE 


0.4 N/kW 
492 days 


4 N/kW 
160 days 
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0.4 N/kW 
518 days 


Pluto 


PLUTO 


4 N/kW 



• 90t spacecraft 


- 50t cargo, 20t power, 
20t propulsion 


• 2MW power 
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1000 AU 


0.4 N/kW 
5.6 years 


4 N/kW 
1.8 years 


90t spacecraft 

] - 50t cargo, 20t power, 
20t propulsion 

2MW power 














PROXIMA CENTAURI 



0.4 N/kW 
122.5 years 



4 N/kW 
29.9 years 


• 90t spacecraft 

- 50t cargo, 20t power, 
20t propulsion 

• 2MW power 
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Original Matthew Jeffries concept from mid 
1960's, rendered by Mark Rademaker 





* 

| 

m 

4 


Rademaker\drex files 


Matthew Jeffries is the artist that created the familiar Star Trek enterprise 
look 
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Warp Field Physics Vision 


Fast Mars 

y 

'. 


Fast Outer 
Solar System 
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Reduced low thrust 
transit times 
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TERRESTRIAL 

PHASE 


Time of Flight 
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Warp Field Interferometer 
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Rapid Interstellar . 


Aggregate negative pressure 
generators in desired topology 
to enable a progression of 
increasingly distant space 
mission destinations 
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Principles of Q-thruster Operation 


Local mass concentrations, say in the form of a 
conventional capacitor with a ceramic dielectric, 
affect vacuum fluctuation density according to 
equation 1 


P 


v local 


= A 


P 


m local 


A 


= JP, 


m local P v 


(i) 


Just as relativistic acceleration (Unruh radiation) can 
change the apparent relative density of the vacuum, 
so too can higher order derivatives according to 
equation 2. 


Sp = 


8p = 


1 




4 nG 

1 


1 ( da\ 1 d 2 a 

+ 


V 


a 


V 


dt J 


( + f AA\ 


4 ttG 


V 




d(j) 

\dt 



a = 


-V(/> 


( 2 ) 


The tools of MagnetoHydroDynamics (MHD) can be 
used to model this modified vacuum fluctuation 
density analogous to how conventional forms of 
electric propulsion model propellant behavior. 















































Gravitational Coupling Constant 



9 Consider the following thought experiment: what would an inertial observer in deep space find if 
the dark energy density were to be integrated over the light horizon of the observable universe, 
~13.7 billion light years? 

9 Starting with the Friedman Equation (and after some manipulation), the following equation can be 
derived that formally captures the results of this thought experiment: 


Vacuum Energy 
Density 



4 7VC t 


2 , 2 

H 



\ 



J 



Equation can be 
rearranged into the 
following form 


Using 9.9xl0' 27 kg/m 3 [2] with t H of 13.7 billion years yields a predicted 
value for the gravitational constant of 6.45x10 11 m 3 /kg*s 2 

A possible physical meaning to this rearranged equation solved for G is 
that gravitation is an emergent phenomenon rather than a fundamental 
force. 

To be specific, the claim could be made that the gravitational coupling 
constant may be a long wavelength consequence (A=ct H ) of dark energy. 



Light horizon 


Dr. Harold “Sonny” White 























Bohr Radius 


* The vacuum perturbation equation just derived can be used to evaluate the state of the 
quantum vacuum in close proximity of the proton at the center of the Hydrogen atom. 

* The first step is to calculate a quasi-classical density for the hydrogen nucleus. The radius of 
the hydrogen atom nucleus is given as /? 0 =1.2x10‘ 15 m {R=R 0 A' /3 where R 0 = 1.2x10‘ 15 m and A is 
the atomic number - these are experimentally determined by electron scattering). 

* The radius can be used with the mass of a proton to calculate a quasi-classical density of the 
hydrogen nucleus: 



• Using p^ 2 / 3 * 9.9x10 -27 kg/m 3 , along with this quasi-classical density p m , the perturbed negative 
pressure state of the quantum vacuum around the hydrogen nucleus is calculated to be: 



The question can be asked how much volume of this perturbed state of the quantum vacuum 
is needed to have the equivalent energy value as the ground state of Hydrogen (13.6eV or 
2.18x1 O' 18 Nm) 




E 


r = 



V 


* The calculated radius is r= 5.29x10 11 m, which is an exact match to the given value for the 
Bohr Radius, a 0 = 5.29x1 O' 11 m. 


Dr. Harold “Sonny” White 














Electron Mass 



Frank Wilczek, Nobel laureate: “We have achieved a beautiful and 
profound understanding of the origin of most of the mass of ordinary 
matter, but not of all of it. The value of the electron mass, in particular, 
remains deeply mysterious...” 

Consider the energy state of the perturbed quantum vacuum field around 
the proton, and set this equal to the kinetic energy of the orbiting electron 
at the ground state. 

4 3 2 1 2 

Pv localC =~ m e V 


We know the speed of the orbiting electron: 

v = ac = c /137 

We can solve for the electron mass, and using the predicted value for 
Pv local °f 3.9x10 -5 kg/m 3 , we get a predicted electron mass of 9.1x10 -31 kg. 


8 


m e = 


m 0 Pv local 0 


f \ 2 

C ' 


V 


137 


/ 
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Magnetic Pressure 


• The first step now is to calculate the magnetic pressure around the Hydrogen 
nucleus. 


The magnetic field as perceived by 
the electron is given by the 
following relationship. The speed of 
the orbiting electron is ac. 


The magnetic pressure is a 


/w 


simple calculation: 

R 




• The quasi-classical plasma pressure of the perturbed quantum vacuum state 
around the Hydrogen nucleus can be calculated by converting the electron velocity 
to temperature using 1/2 m e v 2 =kT , and making the assumption that the virtual 
electron-positron plasma has the same effective temperature as the orbiting 
electron. 

* When the plasma pressure calculation makes use of a 2/3 factor, analogous to the 
predicted dark energy fraction of 2/3 picked up during integration to calculate the 
Gravitational constant, the values are nearly identical: 



Dr. Harold “Sonny” White 














